Mesoscale modeling reveals formation of an epigenetically driven HOXC gene hub.

Gene expression is orchestrated at the structural level by nucleosome positioning, histone tail acetylation, and linker histone (LH) binding. Here, we integrate available data on nucleosome positioning, nucleosome-free regions (NFRs), acetylation islands, and LH binding sites to "fold" in silico the 55-kb HOXC gene cluster and investigate the role of each feature on the gene's folding. The gene cluster spontaneously forms a dynamic connection hub, characterized by hierarchical loops which accommodate multiple contacts simultaneously and decrease the average distance between promoters by ∼100 nm. Contact probability matrices exhibit "stripes" near promoter regions, a feature associated with transcriptional regulation. Interestingly, while LH proteins alone decrease long-range contacts and acetylation alone increases transient contacts, combined LH and acetylation produce long-range contacts. Thus, our work emphasizes how chromatin architecture is coordinated strongly by epigenetic factors and opens the way for nucleosome resolution models incorporating epigenetic modifications to understand and predict gene activity.

Chromatin Fiber Folding Directed by Cooperative Histone Tail Acetylation and Linker Histone Binding.

In eukaryotic chromatin, islands of histone tail acetylation are found near transcription start sites and enhancers, whereas linker histones (LHs) are localized in intergenic regions with wild-type (WT) histone tails. However, the structural mechanisms by which acetylation, in combination with LH binding, modulates chromatin compaction and hence transcription regulation are unknown. To explore the folding propensity by which these features may govern gene expression, we analyze 20 kb fibers that contain regularly spaced acetylation islands of two sizes (2 or 5 kb) with various LH levels by mesoscale modeling. Specifically, we investigate the effect of acetylating each histone tail (H3, H4, H2A, and H2B) individually, in combination (H3 and H4, or all tails), and adding LH to WT regions. We find that fibers with acetylated H4 tails lose local contacts (<1 kb) and fibers with all tails acetylated have decreased long-range contacts in those regions. Tail interaction plots show that this opening of the fiber is driven by the loss of tail-tail interactions in favor of tail-parent core interactions and/or increase in free tails. When adding LH to WT regions, the fibers undergo hierarchical looping, enriching long-range contacts between WT and acetylated domains. For reference, adding LH to the entire fiber results in local condensation and loss of overall long-range contacts. These findings highlight the cooperation between histone tail acetylation and regulatory proteins like LH in directing folding and structural heterogeneity of chromatin fibers. The results advance our understanding of chromatin contact domains, which represent a pivotal part of the cell cycle, diseased states, and differentiation states in eukaryotic cells.

Cohesin Loss Eliminates All Loop Domains.

The human genome folds to create thousands of intervals, called "contact domains," that exhibit enhanced contact frequency within themselves. "Loop domains" form because of tethering between two loci-almost always bound by CTCF and cohesin-lying on the same chromosome. "Compartment domains" form when genomic intervals with similar histone marks co-segregate. Here, we explore the effects of degrading cohesin. All loop domains are eliminated, but neither compartment domains nor histone marks are affected. Loss of loop domains does not lead to widespread ectopic gene activation but does affect a significant minority of active genes. In particular, cohesin loss causes superenhancers to co-localize, forming hundreds of links within and across chromosomes and affecting the regulation of nearby genes. We then restore cohesin and monitor the re-formation of each loop. Although re-formation rates vary greatly, many megabase-sized loops recovered in under an hour, consistent with a model where loop extrusion is rapid.

Kilobase Pair Chromatin Fiber Contacts Promoted by Living-System-Like DNA Linker Length Distributions and Nucleosome Depletion.

Nucleosome placement, or DNA linker length patterns, are believed to yield specific spatial features in chromatin fibers, but details are unknown. Here we examine by mesoscale modeling how kilobase (kb) range contacts and fiber looping depend on linker lengths ranging from 18 to 45 bp, with values modeled after living systems, including nucleosome free regions (NFRs) and gene encoding segments. We also compare artificial constructs with alternating versus randomly distributed linker lengths in the range of 18-72 bp. We show that nonuniform distributions with NFRs enhance flexibility and encourage kb-range contacts. NFRs between neighboring gene segments diminish short-range contacts between flanking nucleosomes, while enhancing kb-range contacts via hierarchical looping. We also demonstrate that variances in linker lengths enhance such contacts. In particular, moderate sized variations in fiber linker lengths (∼27 bp) encourage long-range contacts in randomly distributed linker length fibers. Our work underscores the importance of linker length patterns, alongside bound proteins, in biological regulation. Contacts formed by kb-range chromatin folding are crucial to gene activity. Because we find that special linker length distributions in living systems promote kb contacts, our work suggests ways to manipulate these patterns for regulation of gene activity.

Mesoscale Modeling Reveals Hierarchical Looping of Chromatin Fibers Near Gene Regulatory Elements.

While it is well-recognized that chromatin loops play an important role in gene regulation, structural details regarding higher order chromatin loops are only emerging. Here we present a systematic study of restrained chromatin loops ranging from 25 to 427 nucleosomes (fibers of 5-80 Kb DNA in length), mimicking gene elements studied by 3C contact data. We find that hierarchical looping represents a stable configuration that can effectively bring distant regions of the GATA-4 gene together, satisfying connections reported by 3C experiments. Additionally, we find that restrained chromatin fibers larger than 100 nucleosomes (∼20Kb) form closed plectonemes, whereas fibers shorter than 100 nucleosomes form simple hairpin loops. By studying the dependence of loop structures on internal parameters, we show that loop features are sensitive to linker histone concentration, loop length, divalent ions, and DNA linker length. Specifically, increasing loop length, linker histone concentration, and divalent ion concentration are associated with increased persistence length (or decreased bending), while varying DNA linker length in a manner similar to experimentally observed "nucleosome free regions" (found near transcription start sites) disrupts intertwining and leads to loop opening and increased persistence length in linker histone depleted (-LH) fibers. Chromatin fiber structure sensitivity to these parameters, all of which vary throughout the cell cycle, tissue type, and species, suggests that caution is warranted when using uniform polymer models to fit chromatin conformation capture genome-wide data. Furthermore, the folding geometry we observe near the transcription initiation site of the GATA-4 gene suggests that hierarchical looping provides a structural mechanism for gene inhibition, and offers tunable parameters for design of gene regulation elements.

Design, synthesis, and activity of a series of arylpyrid-3-ylmethanones as type I positive allosteric modulators of α7 nicotinic acetylcholine receptors.

A series of novel arylpyrid-3-ylmethanones (7a-aa) were designed as modulators of α7 nicotinic acetylcholine receptors (nAChRs). The methanones were found to be type I positive allosteric modulators (PAMs) of human α7 nAChRs expressed in Xenopus ooctyes. Structure-activity relationship (SAR) studies resulted in the identification of compound 7v as a potent and efficacious type I PAM with maximum modulation of a nicotine EC5 response of 1200% and EC50 = 0.18 µM. Compound 7z was active in reversing the effect of scopolamine in the novel object recognition (NOR) paradigm with a minimum effective ip dose of 1.0 mg/kg (2.7 µmol/kg). This effect was blocked by the selective α7 nAChR antagonist methyllycaconitine (MLA). These compounds are potent type I positive allosteric modulators of α7 nAChRs that may have therapeutic value in restoring impaired sensory gating and cognitive deficits in schizophrenia and Alzheimer's disease.

Probing sequence-specific DNA flexibility in a-tracts and pyrimidine-purine steps by nuclear magnetic resonance (13)C relaxation and molecular dynamics simulations.

Sequence-specific DNA flexibility plays a key role in a variety of cellular interactions that are critical for gene packaging, expression, and regulation, yet few studies have experimentally explored the sequence dependence of DNA dynamics that occur on biologically relevant time scales. Here, we use nuclear magnetic resonance (NMR) carbon spin relaxation combined with molecular dynamics (MD) simulations to examine the picosecond to nanosecond dynamics in a variety of dinucleotide steps as well as in varying length homopolymeric A(n)·T(n) repeats (A(n)-tracts, where n = 2, 4, or 6) that exhibit unusual structural and mechanical properties. We extend the NMR spin relaxation time scale sensitivity deeper into the nanosecond regime by using glycerol and a longer DNA duplex to slow overall tumbling. Our studies reveal a structurally unique A-tract core (for n > 3) that is uniformly rigid, flanked by junction steps that show increasing sugar flexibility with A-tract length. High sugar mobility is observed at pyrimidine residues at the A-tract junctions, which is encoded at the dinucleotide level (CA, TG, and CG steps) and increases with A-tract length. The MD simulations reproduce many of these trends, particularly the overall rigidity of A-tract base and sugar sites, and suggest that the sugar-backbone dynamics could involve transitions in sugar pucker and phosphate backbone BI ↔ BII equilibria. Our results reinforce an emerging view that sequence-specific DNA flexibility can be imprinted in dynamics occurring deep within the nanosecond time regime that is difficult to characterize experimentally at the atomic level. Such large-amplitude sequence-dependent backbone fluctuations might flag the genome for specific DNA recognition.